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Abstract 

High-energy photons and other particles carrying non-zero orbital angular mo- 
mentum (OAM) emerge as a new tool in high-energy physics. Recently, it was sug- 
gested to generate high-energy photons with non-zero OAM (twisted photons) by 
the Compton backscattering of laser twisted photons on relativistic electron beams. 
Twisted electrons in the intermediate energy range have also been demostrated ex- 
perimentally; twisted protons and other particles can in principle be created in a 
similar way. Collisions of energetic twisted states can offer a new look at particle 
properties and interactions. A theoretical description of twisted particle scattering 
developed previously treated them as pure Bessel states and ran into difficulty when 
describing the OAM of the final twisted particle at non-zero scattering angles. Here 
we develop further this formalism by incorporating two additional important fea- 
tures. First, we treat the initial OAM state as a wave packet of a finite transverse 
size rather than a pure Bessel state. This realistic assumption allows us to resolve 
the existing controversy between two theoretical analyses for non-forward scatter- 
ing. Second, we describe the final twisted particle in terms of the orbital helicity — 
the OAM projection on its average direction of propagation rather than on the fixed 
reaction axis. Using this formalism, we determine to what extent the twisted state 
is transferred from the initial to final OAM particle in a generic scattering kinemat- 
ics. As a particular application, we prove that in the Compton backscattering the 
orbital helicity of the final photon stays close to the OAM projection of the initial 
photon. 



1 Introduction 

1.1 Particles carrying orbital angular momentum 

In high-energy physics we probe the structure of particles and their interactions by bring- 
ing them into collision and detecting the products of their scattering. In general, the more 
control we have on the initial state particles, the more subtle features we can measure. 
For example, sufficiently monochromatic initial beams allow for a direct measurement of 
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the shape of resonances and of interference patterns, while a well defined polarization of 
the initial particles gives access to spin-dependent structure of hadrons. It now seems 
possible that yet another degree of freedom, the orbital angular momentum (OAM) of the 
initial states, can be exploited in the high-energy particle scattering. 

Laser beams carrying non-zero OAM are well known in optics, pQ, for a review see [2], 
numerous applications of light with orbital angular momentum are described in the recent 
book [3]. The light-field in such a beam is described by a non-plane wave solution of the 
Maxwell equations with a helical wave front and an associated integer winding number 
m. Each photon in this light-field, which we call a twisted photon, carries a non-zero value 
of OAM projection onto its propagation axis: L z = mh. An experimental realization jl] 
exists for states with OAM projections as large as m = 200. 

So far, experiments with twisted light were confined mostly to the optical energy range. 
However it was recently noted that Compton backscattering of twisted optical photons 
off an ultra-relativistic electron beam can generate high-energy photons carrying non-zero 
OAM [SI E]. The technology of Compton backscattering is well established [7], thus the 
realization of this idea seems feasible. 

It must be stressed that the possibility to carry OAM is by no mean an exclusive 
property of photons. Other particles can carry orbital angular momentum too. Having 
a non-zero mass or being a fermion does not forbid the existence of phase vortices in 
the transverse plane. Indeed, following the suggestion made in [8], very recently several 
groups have reported successful creation of twisted electrons, first using phase plates \±0\ 
and then with computer-generated holograms [11]. Such electrons carried the energy as 
high as 300 keV and the orbital quantum number up to m ~ 100. It is very conceivable 
that when these electrons are injected into a linear electron accelerator, their energy can 
be boosted into the multi-MeV and even GeV region. Even more, when beams of protons, 
neutrons or other particles with sufficient transverse coherence (and the mere observation 
of neutron diffraction on crystals proves that such coherence is achievable) pass through 
a specially prepared diffractive grating, they can gain an OAM as well. 

One can therefore imagine that with the future progress in this field creation of en- 
ergetic twisted particles of different kind will be possible and can be used in scattering 
experiments. As it was described in [12] , the new degree of freedom that enters this scat- 
tering process might become a new promising tool in nuclear and high-energy physics. 
For example, it can provide access to such features in the structure of hadrons which are 
difficult to probe otherwise. 

1.2 How does OAM change after scattering? 

When a twisted state, be it a photon, an electron or another particle, scatters, elastically 
or inelastically, one can ask how its OAM changes after the scattering. This question 
received little attention so far, mostly because optical photons with OAM are almost 
always assumed to be absorbed rather than scattered. In the case of twisted electrons, 
analyses are limited to semiclassical dynamics in external potentials, see e.g. [5J- The 
fully relativistic, quantum-field-theoretic treatment of this problem, which is absolutely 
necessary for high-energy collisions of twisted particles, has not yet been given. 

In fact, this question is also very crucial for the suggestion of Ej. It was noted that 
for a strictly backward Compton scattering the final energetic photon moving along the 
initial collision axis carries away exactly the same OAM as the initial photon: w! = m. 
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However this conclusion is valid only for a single point of the final phase space. In order 
for this suggestion to become a reliable technique of generation of high-energy twisted 
photons with more or less definite OAM, one must show that m! ~ m holds for small but 
non-zero angles of the final photons 0, at least within the range /3 < m e /E e , where the 
Compton scattering receives its dominant contribution to the total cross section. 

This is the point where a controversy in the literature starts. On one hand, in the 
original paper [5] it was argued, on the basis of an approximate consideration of the non- 
forwards case, that at very small transverse momentum transfer, |P| <C x, the final m' 
stays close to m: 



Here P = p' — p is the momentum transfer to the electron and x ~ (9(1 eV) is the 
conical momentum spread in the initial twisted state. On the other hand, the exact 
non-forward scattering analyzed in [12] for a generic twisted scalar case implies that the 
entire m'-region contributes homogeneously to the cross section at any non-zero transverse 
momentum transfer: 

m! G (-oo, +oo) , P^O. (2) 

These two results are in a clear conflict with each other. 

There are two additional reasons to find these results disturbing. First, if (JT|) holds, 
one can expect that for any reasonable transverse momentum transfer, which is orders of 
magnitude larger than x, the final ml should spread over a very broad range of values 
with almost no correlation with m. If this were true, that would make the experimental 
realization of the suggestion of [3 [6] unfeasible. Second, the result (J2J) is in sharp contrast 
with the conclusion that ml — m for P = 0, which indicates that there is no smooth 
non-forward to forward transition. 

In this paper we resolve all these problems by developing further the formalism of 
twisted particle scattering. First, we allow the initial twisted particles to be more or less 
transversely collimated wave packets, in contrast to the previously analyzed case of pure 
Bessel states. This modification reveals the origin of the discrepancy between (JT|) and 
02]): they correspond to two different limits in the description of non-forward to forward 
transition in twisted state scattering. 

Second, following suggestion of [T2], we describe each OAM state as a twisted state 
with respect to the average propagation axis this very photon rather than using an OAM 
projection on an arbitrarily chosen axis. The OAM projection on the particle averaged 
propagation direction, which we call the orbital helicity, is a more faithful representation 
of the twisted nature of this state. This updated formalism leads us to the remarkable 
conclusion that for small-angle scattering the quantum number ml indeed stays close to 
m even when |P| ^> x. 

In this paper we consider scalar particle scattering with an isotropic matrix element. 
This simplest set up allows us to focus on the universal kinematical features of all scat- 
tering processes, involving a twisted state of a photon, an electron, a proton etc. both 
in the initial and final states. As a particular example, this includes, but is not limited 
to, the Compton scattering of twisted photons. Our generic scalar analysis is related to 
these particular cases just as the scalar theory of diffraction is related to the real diffrac- 
tion of light or of matter waves. In a sense, we investigate what the energy-momentum 
delta-function S(pf — pi) turns into when we pass from plane waves to twisted states. As 
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explained in [Tj2], any non-trivial matrix element will appear as a multiplicative factor in 
front of the resulting expression. 

The paper is organized as follows. In Section [2] we introduce the scalar twisted states, 
remind the reader how scattering of a twisted particle is described, and rederive results 
(DO) and (j2J). In Section [3] we introduce the wave packets into description of twisted states 
and reconcile these results. In Section H] we generalize the formalism to include the orbital 
helicity and finally answer the question of how twisted the final state is. Section [5] contains 
our conclusions. 

Throughout the paper we use the relativistic units K — c — 1. For a 4- vector p = 
(E p ,p), we will separate its 3- vector p into the transverse vector p and the longitudinal 
component p z . Note also that whenever we say forward scattering we actually mean a 
scattering at zero transverse momentum transfer, which might be either strictly forward 
or strictly backward. 



2 Kinematical features of twisted particle scattering 
2.1 Describing twisted states 

Here we briefly summarize the formalism of Bessel-beam twisted states introduced in [5]. 

We first fix a z axis and solve the free wave equation in cylindric coordinates r, (p r , z. 
A solution |x, m) with definite frequency u, longitudinal momentum k z , modulus of the 
transverse momentum |k| = x and a definite z-projection of orbital angular momentum 
m has the form 

|x, m) = e - iut+ik ** ■ ij Km (r) , ^ xm (r) = -=^J m (xr) , (3) 

v 2tt 



where J m (x) is the Bessel function. Such a state is possible for a particle of any mass M; 
the mass appears only in the relation between u, k z and x: u> 2 = k 2 z + x 2 + M 2 . The 
transverse spatial distribution is normalized according to 

/poo 
rf2r C'm'( r )Vw(r) = t^/Vxx 7 J rdrJ rn (xr)J m (x'r) = 5 mm > 5{x - x') . (4) 

A twisted state can be represented as a superposition of plane waves: 

| Xj m) = e - Mk >* J ^a_(k)e jkr , (5) 

where 



. K 

This expansion can be inverted: 

/2tt +oc 

g— iujt+ik z z g*kr 1 -"' 




— i m e- mVk \x,m} , x=|k|. (7) 



x 

m=— oo 



More details about properties of twisted states, their normalization and phase space den- 
sity can be found in [SI [12] • Here we just note that although the wave oscillation am- 
plitude decreases at large radii, the pure Bessel twisted state of finite amplitude is still 
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not localized and not normalizable in the transverse plane. Therefore, the intermediate 
calculations with these Bessel-beam states must be carried out inside a large but finite 
cylindric volume of radius R. 



2.2 Collision of twisted particles 

Let us now re-derive the conflicting results ([1]) and ([2]). 

Consider a 2 — > 2 scattering of plane wave states with initial momenta p and k and 
final momenta p' and k'. The scattering matrix element has the standard form 

t (27r) 4 S^(p + k-p'-k')-M 
S PW (p,k, P} k) = , (8) 

where the amplitude M. is calculated by Feynman rules. The passage from the plane 
wave with momentum k to the twisted state |x, m) can be performed by integrating out 
the plane- wave scattering matrix element over all k with the weight factor a xm (k), as in 
()5]). If we consider elastic scattering between a twisted state and a plane wave, we have 

/d 2 ~k d 2 k.' 
(k)a* Hlm ,(k>)S PW (p,k,p',k'). (9) 

Following [51 [12], we assume in this and the next Section that both the initial and the 
final twisted states here are defined with respect to the common axis z, which is also the 
direction of propagation of the initial plane wave particle with momentum p. 

The integrals in are killed by the transverse delta-function present in 5^ (p + k — 
p' — k'). The presence of a non-trivial amplitude A4, which is a smooth function of the 
momenta, does not influence this integral. Therefore, the key quantity that enters the 
twisted scattering matrix element is the following master integral 

Z w (x,x',P) = j ^^ 'wWa^OO • (27r)V 2 )(k-k'-P) (10) 

where as before P = p' — p is the transverse momentum transfer, and ip and ip' are the 
azimuthal angles of k and k', respectively. It is seen from this expression that the moduli 
of the transverse momenta satisfy the triangle rules (Fig. [1]): 

x<x+|P|, x'<x+|P|, |P|<x + x\ (11) 

After integration over k and |k'|, the master integral can be presented in the form 

Z mm ,(x,x',P) = J\^(m-m 'W-mM 5 (V^l + t 2 + 2t COS ip' - x) d^', (12) 



where 



A = i m '- m e i{m - m '^ p , ip = arctan 6Sm(p , , e = l ~4 (13) 

1 + e cos ip x 
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X 

Figure 1: Triangle with the sides x, x 7 , and |P|. 

and ipp is the azimuthal angle of the momentum transfer P. The differential cross section 
is proportional to |X mm /(x, x', P)| 2 . It is this quantity whose m! distribution and x'- 
dependence determines the OAM properties of the final twisted state. 

In the case of strictly forward scattering, P = 0, one immediately obtains 

X mm /(X, x', 0) = 8{h - X)5mm' ■ (14) 

This result means that in the strictly forward scattering the twisted quantum numbers m 
and x are transferred from the initial to the final particle without any change, [3]. 

The master integral for the non-zero transverse momentum transfer, which was calcu- 
lated in [12], is equal to 



A \/ xx' 

X mm /(x, x', P) = — cos[m'4 - (m - m')5 p ] , (15) 

where A = \>cx! sin 5k is just the area of the triangle with sides x, x 7 , |P| shown in Fig. [U 
and 

V + ^-P 2 \ , /x 2 + P 2 -V 2 



5 k = arccos , 5 P = arccos (16) 

\ zxx' J \ 2x|P| / 

are its angles. One sees that the m'-dependence of the differential cross section comes 
from the oscillating cosine squared: 

|Z mm /(x, x', P)| 2 oc cos 2 [m'5k - (m - m')8 p \ . (17) 

At very large m! this cosine is a strongly oscillating function of the moduli of the momenta 
and it can be replaced by 1/2. Therefore, the contribution of arbitrarily large m! is not 
suppressed. 

As already said above, in an accurate analysis one must keep the cylindric quantization 
volume of radius R large but finite. In this case, the differential cross section receives 
approximately homogeneous contribution from the entire m'-region 



Contribution of each partial wave with a given m' is suppressed by l/m max , so it is the 
summation 



m 



2J |Z w (x, x, P) 



max , _ 



that stays constant in the R — > oo limit. Therefore, in this limit we recover the result 
found in |T2] . 
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On the other hand, one can investigate the small momentum transfer approximation 
inside the master integral. Denoting |P| = ex' and assuming e< 1, one can represent 
expression f JT2|) as 

Z W (x, x', P) « — / e^m-m'W-me^') ^ - x + VeCOS^') dip' . (19) 
27r Jo 

Note that although e < 1, the value of me in the exponential can be large. We then 
perform the formal small-e expansion of the delta-functional under the integral 

5{x' — x' + x'ecosy?') = <5(x' — x) + x'ecosy/ • <5'(x' — x) + o(e) , (20) 

which is valid on a class of sufficiently smooth functions of x or x 7 , and keep only the 
first term. This leads us to the approximate value for the master integral 

X mm /(x, x',P) = A5(x' - x)J m _ m '(me) , (21) 

which clearly shows a smooth transition to the strictly forward/backward case ( THj) . From 
the properties of the Bessel functions, namely that J n (x) is strongly suppressed at x < n, 
one can infer that at small but non-zero transverse momentum transfer P the quantum 
numbers m' stay close to m, and the result ([1]) follows, [5]. 



2.3 The origin of the discrepancy 

A careful inspection shows that discrepancy between the results ([1]) and ([2]) is linked to 
the transverse spatial extent of the incoming twisted state. 

First of all, as it was already discussed in [12], the discontinuous non- forward to 
forward transition is a consequence of the infinite transverse size of the pure Bessel-beam 
state. If the radius R of the quantization volume is kept large but finite, then the smooth 
transition is restored in an extremely narrow region of |P| ~ 1/R. This suggests that if 
one replaces the pure twisted initial state |x, m) with a wave packet 

POO 

|z) = / dx/(x)|x,m>, (22) 
Jo 

with a narrow weight function /(x) peaked at x = x and having a width a <C x , then 
the non-forward to forward transition is expected to be smooth even for infinite R and to 
take place within the transverse momentum transfer region |P| ~ a. 

On the other hand, the derivation of ([T|) just reproduced involves manipulation of the 
delta-functional, which is valid only if a convolution with a sufficiently smooth function of 
x or x' is assumed. The pure Bessel states lack such a convolution, therefore this result 
cannot be expected to hold for the pure Bessel states. However, for a sufficiently compact 
wave packet this conclusion can be valid. 

These observations necessitate a careful re-analysis of the master integral for a situa- 
tion when the initial state is described by a wave packet (I2"2l rather than a pure Bessel 
state. 
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3 Scattering of a twisted wave packet 



3.1 Qualitative features 

Let us start by discussing qualitative features of the coordinate space wave function of 
the wave packet ( l22l . When quantitative estimates are needed, we will use the gaussian 
approximation for the weight function /(x): 

\2" 



/(x) = iVexp 



yX~ Xpj 

2a 2 



(23) 



with a <^ kq and the normalization coefficient N fixed by 

POO 

/ |/(x)| 2 dx=l. (24) 
Jo 

The transverse coordinate wave function ip m ( r ) of this wave packet is defined by 

dxf{K)il) H7n {r) = —= / \fxdxJ m (xr)f(x) = —=tfj m (r) (25) 



io v27r J V27T 

and is normalized to unity: 

POO 

rdr\^ m (r)\ 2 = 1- (26) 

In Appendix we study some properties of this averaged wave function. We show there 
that if m is not too large, m <C x 2 /cr 2 , the wave function ip m (r) exhibits radial oscillations 
characteristic of the Bessel function J m (xor) until r becomes larger than the coherence 
radius r c = 1/a. Beyond this radius, the averaged wave function is strongly suppressed. 
This behavior is well seen in Fig. [2] where we plotted ip^{r) for a = x /5 and a = x /20 
for x = 1 (in arbitrary units) and compared it with the pure Bessel state y^J 5 (xor). 
However, if m ^> x 2 /cr 2 , there is no room left for the radial oscillations, and the wave 
function is strongly peaked at r = m/x . 

The effective regularization of the radial wave function by the coherence radius r c = 
1/a plays an important role in the master integral and its m'-dependence. The master 
integral (fTTT) for the pure twisted states can be also represented as a triple-Bessel integral, 
see P2]: 

POO 

2w(x, x',P) = A\[kk' \ rdr J m (xr) J m /(xV) J m _ m /(|P|r) . (27) 

Jo 

If m! is extremely large, the main contribution comes from the large r-region, r > 
max(m'/x', m'/|P|). A pure Bessel function for the initial state, J m (xr), is not suffi- 
ciently suppressed at such large r. However, if the initial state is a wave packet, then 
the exponential suppression is at work beyond r c , which effectively limits the values of 
\m — m'\. 



3.2 Averaged master integral 

Since the master integral is a linear functional of the initial wave function, the averaged 
master integral can be represented as 

POO 

Xw(x', P)= / <ixZ mm /(x, x', P)/(x) (28) 
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Figure 2: Radial wave functions for m = 5 and x = 1: the pure Bessel state ^JJcq Jbi. >c o r ) 
(dash-dotted line), the function ip m (r) for a = 1/20 (solid line), the function ip m {r) for 
a = l/5 (dashed line). 
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Figure 3: The weight function as a function of x in the narrow peak (left) and broad peak 
(right) regimes. 



or, using Eq. (1121) , as 



-^mm' )P) 



2Wo 



2tt 



3 i[(m— m/)ip'— mtf>] 



f (xVi + e 2 + 2ecos< / /) 

/ 2 M7^ d ^'- 

(l + e 2 + 2ecosy/) 7 



(29) 



In this Section we aim at resolving the discrepancy between the results ([I]) and ([2]). 
Therefore, it is sufficient to consider non-zero but very small momentum transfer, |P| = 
ex' with e < 1. In this case x changes in the small interval from x 7 — |P| to x 7 + |P| and 
the averaged master integral can be approximated as 



/ ,P) 



A 
2^ 



2tt 



j[(m—m')ip' —em sin if>'\ 



f (x' + |P| cosy?') dif' . 



(30) 



For the further analysis we need to distinguish two cases depending on the shape of 
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f{x l + |P| cosy?') as a function of x = cosy?' on the interval [—1,1]: 
/(x' + x|P|) = iVexp 



(x - x ) 2 



,T ° = ipi ' X = W\' ( ' 



2AI 

If it is strongly peaked as in Fig. [3] left, we are dealing with the narrow-peak situation. 
It corresponds to |xo| < 1 and A x <C 1: 

narrow peak: a <C |P| (32) 

In this case the relevant integrals can be taken approximately using Laplace method. The 
opposite case, A x ^> 1, corresponds to the broad-peak situation shown in Fig. El right: 

broad peak: |P| <C a <C x . (33) 

In this case the integrals can be approximately calculated by using the Taylor expansion 
of the weight function. 

Let us derive the m! distribution in the broad and narrow peak situations. In the 
broad peak the first approximation, one can put / (x' + |P| cosy/) ~ / (x 7 ) and 

immediately obtain the result 

Z mm ,(x', P) « Af (x') J m _ m ,(em) . (34) 

Using the properties of the Bessel functions, one can deduce that m' should always stay 
close to m, otherwise the contribution is suppressed. For small values of m, when me <C 1, 
the only significant contribution comes from m! = m. For large m, when me 3> 1, the 
m'-distribution is spread over several values around m' ~ m. 

For the narrow peak case, \m — m'\ 3> me still holds almost in the entire m! region. 
Then, we can read off ( 130]) that the averaged master integral effectively extracts the 
(m — m')-th Fourier harmonic of the weight function. Since the weight function has a 
peak with a width A x , we conclude that 

1 IPI 

\ m -'> ri \% — = — . (35) 
A x a 

This means that the final m! can strongly differ from and be much larger than the initial 
m if a is sufficiently small. Considering additionally the region \m — m'\ < me does not 
change this conclusion. 



3.3 Different limits 

The key quantity which we discuss in this Section is the m'-distribution of the scattering 
matrix element at small P. Our analysis reveals that the answer to this question depends 
in fact on a subtle interplay between two different limits: a — > and |P| — > 0. The 
apparently conflicting results ([T|) and ([2]) simply correspond to two different choices of 
which limit is taken first. 

If the initial wave packet is fixed (i.e. a is kept constant) and |P| — > 0, then we are in 
the broad-peak regime, and the result (JT|) follows. In this way we study the non-forward 
to forward transition for a wave packet of fixed size. This implicit assumption, although 
not mentioned in [5], is essential when deriving this result. 

On the contrary, if |P| is kept constant but o decreases, we enter the narrow-peak 
regime, and fl35l) holds. In the a — > limit we recover the result (T5]). This limit corresponds 
to the point- by-point analysis of the small-|P| behaviour in the scattering of true Bessel 
states, studied in |l2j . This resolves the discrepancy mentioned in the introduction. 
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4 Orbital helicity of the final twisted state 



In the previous Section we reconciled results ([I]) and ([2]) which was a technical rather 
than physical problem. However our calculations did not answer the physically important 
question: what are the true OAM properties of the final twisted state with respect to its 
own propagation direction? 

The results of the previous Section imply that if the momentum transfer |P| ^> x, 
a very broad m'-region contributes to the differential cross section. However, this broad 
region can easily be an artefact of using the same z axis to describe both the initial and 
the final twisted states. Indeed, ([7]) shows that a non-forward plane wave, when expanded 
in the basis of twisted states with respect to axis z, involves all values of m from minus to 
plus infinity, despite the fact that it actually carries no OAM at all. As it was suggested 
in [12] , for the full resolution of this problem one needs to introduce the concept of orbital 
helicity: the projection of orbital angular momentum not on the reaction axis z but on 
the axis of the average propagation direction of the final particle. This is what we do in 
the present Section. 

4.1 Orbital helicity distribution: pure Bessel states 

Let us consider again the scattering process 

\^m) H +\PW(p)} -> |x',m%, + \PW{f)) (36) 

with fixed x, (kn) = k z , m and fixed 3- momentum transfer P = p> — p. The twisted 
states in ( 136]) are defined with respect to two different axes. The incoming state is written 
with respect to the direction n which coincides with the axis z. The final twisted state is 
defined with respect to another direction n' ^ n to be defined below. The average value 
of the 3-momentum in the initial and final states are 

(k) = (x,m\k\x,m) = k z n , (Ar) = (x',m'\k'\x',m') = k' z ,n' . (37) 

Note that k' z , means here the component of the vector k! along the axis z' directed along 
n' . 

A priori, any choice of the axis n' is allowed. We find it convenient to use the following 
prescription: we draw n' along the direction of the well-defined 3-momentum q = {k)—P = 
qn' . Note that although q and (Ar) are parallel, they are not equal, q ^ k' z ,, because the 
averaged momenta are not supposed to obey any conservation law 

(k)^(k')+P. (38) 

In this way the quantity m! becomes the orbital helicity rather than just OAM projection 
on an arbitrary axis. 

We now aim at evaluating the master integral for this scattering. When calculating 
the scattering matrix element for the process ( 136]) . we have a very familiar expression for 
the master integral: 

Zmm' = 27T 8{k z — k z — P z ) • X mm i(^K, K , P) 

= i m '" m v / ^x 7 / ^Ve i(mp " mV '' • o~ {3) (k -k'-P). (39) 
Jo 
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Notice two important differences with respect to the "coaxial" master integral (llip . First, 
the azimuthal angles ip and if' lie in two different planes which are orthogonal to n and n', 
respectively. Second, the master integral contains 3-dimensional delta-function instead of 
just two-dimensional: since the integration is not limited to a single transverse plane, the 
longitudinal delta-function is, naturally, included in the integral. 

The integral fl39|) can be easily done in the usual cartesian coordinate frame (x,y,z). 
This frame is fixed in the following way (see Fig. H]). The directions n and n' define the 
axis z and the (x, z) plane. The axis y is orthogonal to the plane containing n, n', P. In 
this coordinate frame the momenta k and P are 




P 



X 



P= I I , (40) 
P z 

with P x < 0. The direction of n' is (sin /3, 0, cos /3), which is defined by 

/ -P x \ ( sin/3 \ 
q=(k)-P= \=q , (41) 







( sin/3 \ 








\ h - P z ) 




I cos/3 / 



where q 2 = P x + (k z — P z ) 2 . Then, the momentum k! has components 

(x'c^cy + k' z ,sp \ 
x's^ , (42) 

-Vsgcy + k' z ,Cp J 

where we introduced the obvious short-hand notation for sines and cosines. 

The three-dimensional delta-function in (1391) expresses the conservation of the three- 
momentum at the level of plane waves: k = k! + P. Writing it explicitly, we obtain 

S($(k-k>-p) = 5(xc^-x>cpc^-(k' z ,-q)sp) 

x 8{xs v - ds^i) ■ S^'spCpi - {k' z , - q)cp) . (43) 

In order for the integral to be non-zero, we require that 

Aa > \k' z , - q\cp . (44) 
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Then the z-delta-function can be used to kill the ip' integration. Let us for the moment 
drop the exp(irmp — im'<p') factor from the master integral (139]) . Then we have 



d(pd(p'5 {3 \k- k' -P) 
y5[x' 2 + (k' z ,-q) 2 -x 2 ] J d^5^ 



K> - q 

XC, 



<P 



a/ x /2 — x 2 cos 2 f3 



5[x' 2 + (^-g) 2 -x 2 ] . (45) 



Here we implicitly assumed that the integration over ip goes not around the full 2ir domain, 
but only over the semicircle where sirup has the same sign as simp'. The extra delta- 
function in (j45p replaces S(k z — k' z — P z ) in the "coaxial" case. Here it leads to important 
conclusions that 

x' < x, \k' z , — q\ < x. (46) 

So, the values of the momentum transfer P (both P x and P z ) can be large, \P x ,z\ ^ x j 
but they must be accompanied by a correspondingly large value of k' z ,. 
It is convenient to introduce the "angle" £ by 

x' = xcos£, k' z , - q = xsin£, £ G [— 7r/2, tt/2] . (47) 

In this notation condition (jHJ) becomes 

sin(/3 - |e|) > ^ /3>|e|- (48) 

which also implies x 2 > x' 2 > x 2 cos 2 (5. The pair of values of <p and <p', which are set by 
the delta-functions, can be written as 

sin£ , tan£ 

cos(p = ^—, cos(p=- -. 49 

sin p tan p 

Let us now recall that the master integral (13 9 p contains the exponential factor exp(im(p— 
im'ip'). Since the integral just taken receives its contribution only from two (ip, ip') points, 
this extra factor is simply an overall multiplier computed at each of these two points. 
Effectively, it corresponds to the replacement 2 — > 2cos(rrup — rn'ip'). Thus, the final 
result for the master integral is 

X mm> = ,-W^7 ■ 2 c f™V - ™V) 5[h - + {K) _ q f _ x 2 ] , (50) 

a/ K a — K l COS^ p 

where the values of ip and ip' are given by (j490 . This expression shows the final orbital 
helicity m! distribution. Similarly to the "coaxial" case, X mm / is an oscillatory function of 
m! . Clearly, |X mm /(x, x')| 2 is not suppressed even at extremely large m'. However, this is 
an artefact of taking pure Bessel beams for the initial and final twisted states. In reality, 
a pure Bessel state is as unphysical as a plane wave, because its radial coordinate wave 
function effectively extends to infinity and is not normalizable. Even if m is small, all 
values of m' up to infinity contribute to the cross section. 
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4.2 Orbital helicity distribution: wave packets 

The above conclusion is expected to change if the initial and final twisted states are as- 
sumed to be a wave packets rather than pure Bessel states. As before, we represent the 
initial and final twisted states similarly to ( 125]) : initial twisted state as a superposition of 
\x, m) states with equal value of m and different values of x (distributed by the weight 
function f(x) around x in a region with width a) and final twisted state as a superpo- 
sition of |x',m') states with equal value of m' and different values of x 1 (distributed by 
the weight function g(x') around Xq in a region with width a'). 
Then we calculate the averaged master integral 



-yav 
mm' 



J dxf(x) J dx'g(x')l mm/ (x,x') . (51) 



Since the initial laser photon is assumed to be monocromatic, the final photon has a fixed 
energy, oj' = E + oj — E', therefore, the longitudinal momenta 

k z = Vu 2 - x 2 , k' z , = vV) 2 - (x') 2 (52) 

change when x and x' are varied. 

Let us first use the x' integration to eliminate the delta-function. Manipulation with 
delta-function gives: 



q 2 + u' 2 -x 2 ( , 2D\ 



S \x a + [k' z , - qf - X 2 ] = 6 (x 2 - ^ -q 2 + 2,Vw° - *") = 9 + * 6 ( 

where 

D = - ^2 (x 2 co' 2 + cu' 2 q 3 + q 2 x 2 ) - x 4 - co' 4 - q 4 . 
After the x 1 integration we obtain 

ir„, = / f-M ■ 2 c ° s ^- 2 mV 2 \ i!±" ( 53) 

J \ x a/x /2 - x 2 cos 2 (3 4g 

with x 7 = 2D/g. 

The resulting expression for X™ m , can be easily analyzed in the case of small-angle 
scattering, ^ < 1, which implies |£| < /3 -C 1 and <£> ~ y/. As the result, the cosine 
becomes simply 

cos[(m - m!) arccos(£//3)] _ r m /_ m (£/ft) 

VCOS 2 £ - COS 2 J) ~ y//3 2 - £ 2 ' 

where T n (x) is the Chebyshev's polynomial of the first kind: 

cos(n • arccos(a:)) = T n (x) . 
The integral is then proportional to 

dxf{x)g{x') Tm ' ~ m ^/^ = f ^ F (0 Tm, ' m(e//3 ) , (55) 

where F(£) includes both weight functions and the Jacobian of the x — > £ change of 
variables. 



14 



This Chebyshev polynomial makes (m/ — m)/4 oscillations on the interval £//3 from 

— 1 to 1. Smearing of x due to the weight function leads to smearing of by the 
amount of A£/(8/3). Therefore, if |m' — m| ^> 8/3/A£, the oscillations strongly suppress 
the contribution. We conclude that only ml such that 

I™ ~ m l £ ^7 ( 56 ) 

effectively contribute to the cross section. 
Let us estimate A£ using expression 



t « (57) 
x 

obtained from Eq. (f471) . From here we have 

. , A A;', -Aq k',-q. Aq Ax a 

Af« — 5 - z 9 Ax^ -» (58) 

x x x u; a; 

since the following estimations are valid from (152|) at w < w' and x 7 « x: 

xAx . ,. VAx' . , . . 

AA; 2 w , AA^, w — < Afc 2 , (59) 

and therefore, 

Aq = Ay/(k z - P z f + Pi = ih - P ^ z)Akz » -- Ax. (60) 
As a result, we obtain the important estimate: 



\m'-m\<8(3-. (61) 

Note that the right hand side of this inequality can be small for small enough angles /3. 

A crucial observation is that the result (l6li) does not require the transverse momentum 
transfer |P| to be small (smaller than kq or even a as it was the case in the previous 
Section). The only assumption for this simple analysis was that the scattering angle was 
much smaller than one. We remind once again that for the Compton scattering, where 
the region of small scattering angles (3 ~ m e /E e gives the dominant contribution to the 
total cross section, this approximation is fully valid. 

Therefore, we proved that if the initial and final twisted states are the wave packets, 
then final orbital helicity m! stays close to m and the final x' stays close to x during 
small-angle scattering. 



5 Conclusions 

Photons carrying non-zero orbital angular momentum are now routinely produced in the 
low-energy domain. It was recently argued in [51 [6] that Compton backscattering of 
photons with OAM can help create high-energy photons carrying large values of OAM. 
Intermediate energy electrons with large OAM have been recently created experimentally, 
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and the same technique can in principle be used to give twist to other particles, such as 
protons. Thus, energetic particles with OAM and their scattering can become a new 
promising tool in experimental subatomic physics. 

In these circumstances, it is essential to understand quantum-field-theoretically how 
the OAM of a particles changes after its scattering. In particular, the vital question 
about the OAM of the final twisted photon scattered at non-zero angle, which arises 
in the context of Compton backscattering, has not been satisfactorily answered so far. 
Two previous analyzes, [5] and [12], led to conflicting results. In order to resolve this 
controversy and to answer the question, one needs to generalize the formalism and to 
allow the initial twisted state to be in the form of a wave packet rather than a pure 
Bessel state. This generalization of the formalism is done in the present paper. Using 
it, we completely resolve the existing controversy by showing that the two apparently 
conflicting results emerge in two different limits of the same expression. 

Another key part of this work is the concept of orbital helicity of a twisted particle: 
that is, the OAM projection on its averaged propagation direction rather than on the 
reaction axis. With this more physically appealing definition we derived the orbital helicity 
distribution of a generic scalar scattering at any angle. Our results confirm the intuitive 
expectation that the orbital helicity is approximately conserved at small-angle scattering. 
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A Wave packets of Bessel states 

Here, for the sake of completeness, we study the radial wave function (125 p of the Bessel- 
beam wave packet in some detail. 

Let us start with its expected qualitative behaviour if m is not too large. As r grows, 
the radial oscillations of different twisted states in the wave packed remain almost in phase 
up to the coherence radius r c = 1/a, and at r > r c destructive interference strongly sup- 
presses the wave function. However, at extremely large radii the true large-r asymptotics 
determined by the small- x behavior of the weight function sets in. 

The calculations below show that this qualitative picture changes when m > x 2 /cr 2 . 
Depending on whether m is larger or smaller than x 2 /<7 2 , we will speak of the small-m 
and large-m regions. 

To study the shape of ip m (r), we break it into three r-regions: the central hole, x r < 
m, the intermediate region, x r ^> m, and the true large-r asymptotics, which sets in at 
r > r* to be determined below. In the first two regions, the approximate forms of the 
Bessel functions allow us to represent ip m (r) as 



q 2 (m + 1/2) 



2" 



exp 



at 



x r < m 



ip m (r) w A^v / 27ro- v /xbJ m (x r) x 




(62) 



exp 



at 



Xor 3> m . 
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These shapes come from the region under the peak in the weight function /(x) which is 
assumed to be its dominating feature. 

The very far asymptotics of ip m {r) comes from the small- x behaviour of /(x). For 
example, if /(0) is finite, then /(0) can be taken out of the integral, and we are left with 
an integral of the form of 



The large y behavior of this intregral is 



dx\/xJ n 



x 



(63) 



£(! + !)_ /2 

r(i + f) Vvr 



cos I y 



7T 



-m 



7T 



4 



+ 0(y 



The first term is very close to y/m at all m > 1, while the second term is a small amplitude 
oscillatory function. Thus, a reasonable estimate for the large-r asymptotics of the wave 
function is 



ip m {r) « AT exp 



Xn 



2a 2 



-3/2 



(64) 



where the exponentially suppressed factor is simply /(0). Note that the exact r-power 



-3/2 : 



is a 



is driven by the /(x) behaviour in the x — )■ limit; the asymptotics ip m ( r ) ^ r 
consequence of the finite /(0). 

By comparing f|62|) and (|64p . one finds that the true large-r asymptotics sets in beyond 



r* 



x ^ l 

» r c = - 



a 



(65) 



The intermediate region exists if r* ^> m/x . Therefore, the qualitative picture described 
above holds if m > Xq/ct 2 , that is, in the small-m region. 

For very large m, there is no room for the intermediate r-region. The radius at which 
the Bessel function J m is supposed to have its first peak is already so large, that destructive 
interference takes place. Therefore, the wave function exhibits no radial oscillations and 
has the shape of a single peak at r = m/x . Its shape can be estimated to be 



exp 



2r 2 a 2 



x, 



3/2 



m 



exp 



(r — m/xo) 1 
2A 2 



(66) 



with A = ma/xQ (so that A/r = cr/x ). 
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